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"Is there a q-analogue of that?" — Dennis Stanton. 

Abstract. We prove a Matrix- Tree Theorem enumerating the spanning trees of a cell complex 
in terms of the eigenvalues of its cellular Laplacian operators, generalizing a previous result for 
simplicial complexes. As an application, we obtain explicit formulas for spanning tree enumera- 
tors and Laplacian eigenvalues of cubes; the latter are integers. We prove a weighted version of 
the eigenvalue formula, providing evidence for a conjecture on weighted enumeration of cubical 
spanning trees. We introduce a cubical analogue of shiftedness, and obtain a recursive formula 
for the Laplacian eigenvalues of shifted cubical complexes, in particular, these eigenvalues are also 
integers. Finally, we recover Adin's enumeration of spanning trees of a complete colorful simpli- 
cial complex from the cellular Matrix- Tree Theorem together with a result of Kook, Reiner and 
Stanton. 



1. Introduction 

1.1. Cellular spanning trees. In [8], the authors initiated the study of simplicial spanning trees: 
subcomplexes of a simplicial complex that behave much like the spanning trees of a graph. The 
central result of [5] is a generalization of the Matrix- Tree Theorem, enumerating simplicial spanning 
trees in terms of eigenvalues of combinatorial Laplacians. In this paper, we extend our field of 
inquiry to the setting of arbitrary cell complexes and their Laplacians. 

Let X be a d-dimensional cell complex; we write Xj for the set of i-cells of X . A cellular spanning 
tree, or just a spanning tree, of X is a d-dimcnsional subcomplex Y C X with the same (d — 1)- 
skeleton, whose (reduced) homology satisfies the two conditions Hd{Y',^) = and \Hd-i{Y;li)\ < 
oo. These two conditions imply that \Yd\ — \Xa\ — $d{X) + f3d-i(X), where $d and (3d-i denote 
(reduced) Betti numbers. In fact, any two of these three conditions together imply the third. In 
the case d = 1, a cellular spanning tree is just a spanning tree of X in the usual graph-theoretic 
sense; the first two conditions above say respectively that Y is acyclic and connected, and the third 
condition says that Y has one fewer edge than the number of vertices in X . 

Let Ci denote the i th cellular chain group of X with coefficients in Z, and let d t : Ci — > Cj_i and 
d* : Ci-i — > Ci be the cellular boundary and coboundary maps (where we have identified cochains 
with chains via the natural inner product). The i th up-down, down-up and total combinatorial 
Laplacians are respectively 

^ = ^•+1^1. ^ = 0*8,, L\ ot = L" d + Lf n , 
which may be viewed either as endomorphisms on Ci, or as square symmetric matrices, as convenient. 
We are interested in the spectra of these Laplacians, that is, their multisets of eigenvalues, which 
we denote by s" d (X), sf u (X), and s- ot (X) respectively. Combinatorial Laplacians seem to have 
first appeared in the work of Eckmann [12] on finite-dimensional Hodge theory, in which the i th 
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homology group of a chain complex is identified with ker(L,) via the direct sum decomposition 
Ci — imd i+1 ®ker Li ffiim d* . As the name suggests, the combinatorial Laplacian is a discrete version 
of the Laplacian on differential forms for a Riemannian manifold; Dodziuk and Patodi [5] proved 
that for suitably nice triangulations, the eigenvalues of the combinatorial and analytic Laplacians 
converge to each other. 

For < k < d, let Tk(X) denote the set of all spanning fc-trees of A (that is, the spanning trees 
of the fc-skeleton of X). Let 

r k (X)= \H k -i{Y--L)\\ 
YeT k (x) 

TT k {X) = Yl A for fc > 1. 

We set 7To = |Xo|, the number of vertices of X, for reasons that will become clear later. Note that 
T\(X) is just the number of spanning trees of the 1-skeleton of X . Bolker [4] was the first to observe 
that enumeration of higher-dimensional trees requires some consideration of torsion; the specific 
summand \Hk-i(Y; Z)| 2 , first noticed by Kalai [16], arises from an application of the Binet-Cauchy 
theorem. The precise relationship between the families of invariants {irk(X)} and {rfc(A)} is as 
follows. 

Theorem 12.81 [Cellular Matrix- Tree Theorem]. Let d > 1, and let X d be a cell complex such 
that Hi(X;Q) = for all i < d. Fix a spanning (d — l)-tree Y C X, let U = and let Ljj be 

the matrix obtained from L^ 1 (X) by deleting the rows and columns corresponding to U. Then: 

M X) = T ^ )Trf - lW and 
\H d _ 2 (Y;Z)\ 2 

These formulas are identical to those of the Simplicial Matrix- Tree Theorem [5J Theorem 1.3], 
except that they apply to the wider class of cell complexes. Note that the classical Matrix- Tree 
Theorem for graphs is just the special case d = 1. Our proof, in Section [21 is similar to that in 
[5], and ultimately stems from the work of Kalai [16] . It is not hard to generalize the result to a 
weighted version, Theorem l2.12l 

Of particular interest to us is the problem of which complexes are Laplacian integral, that is, 
all of whose Laplacians have integer eigenvalues. This question seems to have been first raised by 
Friedman |13] . While no general characterization of Laplacian integrality is known, many classes of 
simplicial complexes are known to have combinatorially meaningful Laplacian integer spectra. These 
include matroid complexes [19] , shifted complexes [7], matching complexes [TU], and chessboard 
complexes [14] . 

In Section [5] of the paper, we develop this machinery more carefully and prove the Cellular 
Matrix- Tree Theorem and its weighted analogue. 

1.2. Cubical complexes. In Section [3] we study Laplacian integrality for cell complexes that are 
cubical rather than simplicial. Specifically, the n-dimcnsional unit cube Q n C R" can be regarded as 
a cell complex whose faces are indexed by the ordered n-tuples / — (fx, . . . , /„), where fi € {0, 1, ©}, 
with / C g iff fi = g t whenever gi ^ ©; here dim / = = ©}. Intuitively, we think of © as 

the open interval (0, 1). A proper cubical complex is then an order ideal in this face poset. (This is a 
more restricted class than the more usual definition of "cubical complex" , meaning any cell complex 
in which the closure of each cell is a cube.) 

For a proper cubical complex X C Q n , the prism over X is the cell complex PA :— X x {0, 1, ©} C 
Qn+i- The cellular Laplacian spectra of A determine those of PA; the precise formulas are given 
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in Theorem 13.31 Since the cube Q n itself is an iterated prism, it follows from these results that Q r , 
and all its skeletons are Laplacian integral, with spectra given by the formulas 



E 



Aest ot (Qr») 



r*= "^(l + rT*. (1) 



e ^ e (:)C;V ( 2 ) 



(Theorem l3.4j) . where = means "equal up to constant coefficient" (i.e., up to the multiplicity of zero 
as an eigenvalue) . Applying the Cellular Matrix- Tree Theorem to formula @ gives 

n 

n(Qn)= II m mzX) 
j=k+i 

(Corollary I3.5j) . This generalizes the well-known count rij=2(^i) °f spanning trees of the n-cube 
graph (see, e.g., [24j Example 5.6.10]). 

1.3. Weighted spanning tree enumerators. In Section IU we discuss weighted spanning tree 
enumerators of skeletons of cubes. 

To each face / = (/i, . . . , /„) 6 Q ni associate the monomial 

£/ = n « n Xi n v*> 

»:/i=® »:/i=0 i:/*=l 

in commuting indeterminates. For a cubical subcomplex X C Q„, we consider the weighted spanning 
tree enumerator 

ux)= e no 

YeT k (X)f£Y k 

In principle, one can calculate fk(X) by replacing the Laplacians of X with their weighted versions 
and applying Theorem 12.121 In practice, the difficulty is that the eigenvalues of the weighted 
Laplacians of cubical complexes — even for such basic cases as skeletons of cubes — are in general 
not polynomial. Furthermore, the corresponding weighted cellular boundary maps d do not satisfy 
d 2 = 0, and so do not give the structure of an algebraic chain complex. 

In order to overcome this difficulty we replace the natural "combinatorial" weighting £y with a 
closely related "algebraic" weighting by Laurent monomials, whose boundary maps do satisfy d 2 = 0. 
The algebraically weighted Laplacian spectra of Q n is given explicitly as follows (Theorem 14. 2\i : 



Aes l fc ot (Q™) j=i x ' 

where itj = qf jx 2 + qf /yf and ej denotes the j th elementary symmetric function. 

In the case of shifted simplicial complexes, a weighted version of the Simplicial Matrix- Tree 
Theorem can be used to translate the knowledge of the algebraically weighted Laplacian spectrum 
into an exact combinatorial formula for the weighted tree enumerators tj.(X) (see [8]). We had 
hoped that this approach would succeed in the case of skeletons of cubes, but the exact form of the 
relationship between the two weightings has so far eluded us. On the other hand, there is strong 
evidence that the combinatorial spanning tree enumerator Tk(Q n ) is given by the following formula 
(Conjecture I4.3P : 
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v- 1 fn-l\fi-l\ (\ A \~^ 

k- 1 , 



r fc (g„) = (gi---g„)- fe - 1 J] 



AC[n] 
|A|>fc+l 



1.4. Shifted cubical complexes. In Section^ we introduce an analogue of shiftedness for cubical 
complexes. 

Shifted complexes are a well-known and important family of simplicial complexes; see, e.g., [3"lll7] 
as general references. They may be defined as certain iterated near-cones, or as order ideals with 
respect to a natural partial ordering of integer sets. The Laplacian spectrum of a shifted simplicial 
complex A can be expressed recursively in terms of the spectra of its link and deletion [JJ Lemma 5.3] , 
and explicitly as the transpose of its vertex- facet degree sequence [Jj Theorem 1.1]. In particular, 
the eigenvalues are integers. 

Our definition of a cubical analogue of shiftedness seeks to extend all these properties of shifted 
simplicial complexes to the cubical setting. First, shifted sub-complexes of Q n are in bijection with 
the shifted simplicial complexes on n vertices, via the "mirroring" operation as described in Babson, 
Billera, and Chan [2]; roughly speaking, vertices of a simplicial complex correspond to directions of 
a cubical complex. Like their simplicial cousins, shifted cubical complexes are a particular type of 
iterated near-prism. Moreover, their eigenvalues may be described recursively in terms of the cubical 
link and deletion. An immediate corollary of this recurrence is that shifted cubical complexes are 
Laplacian integral. We present a number of natural questions for further research on shifted cubical 
complexes. 

1.5. Duality. In Scction[6l we study the connections between the Laplacian spectra of a cell complex 
and its dual complex. 

Two cell complexes X, Y are dual if there is an inclusion-reversing bijection /•->■/* from the cells 
of X to the cells of Y, with dim / + dim f*=d for all /, such that the boundary maps of X equal 
the coboundary maps of Y. In this case, we have an immediate duality on the Laplacian spectra: 
sl ot (X) = s^-i(Y) for all i. Moreover, for subsets T C X. t and U = {/* : / G X. t \T}, the subcomplex 
Xt = TUX((_i) forms a spanning i-tree of X if and only if the subcomplex Y\j = [/UYrd-i-i) forms 
a spanning [d — i)-tree of Y (Proposition 16. 1 \ H From a matroid point of view, this is essentially the 
statement that the matroid represented by the columns of d x ■ is dual to the matroid represented by 
the columns of dy^-; the bases of these matroids are cellular spanning trees. (If X is a simplicial 
complex, this matroid is known as a simplicial matroid; see, e.g., [5].) 

The cellular dual of the n-cube is the n-dimensional cross-polytope. This is an instance of a 
complete colorful complex, whose spanning tree enumerators were calculated by Adin [1 . In addition, 
complete colorful complexes are matroid complexes (corresponding to products of rank-1 matroids). 
Kook, Reiner and Stanton determined the Laplacian spectra of matroid complexes in 19 . We 
explore the connection between all these results. In particular, we explain how to derive Adin's 
formula from the Kook-Reiner-Stanton formula together with the Cellular Matrix- Tree Theorem, 
and give an alternate calculation of the (unweighted) spectra of cubes (Theorem 13.41) using duality 
together with Adin's formula. 

We thank Gil Kalai for suggesting cubical complexes as fertile ground for exploration; Eran Nevo 
for solving Problem 15. 8( and two anonymous referees for their careful reading and many valuable 
suggestions. 



As pointed out by a referee, Yjj = {/* : / ^ Xt}, a relationship reminiscent of Alexander duality for simplicial 
complexes. 
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2. Cellular spanning trees 

2.1. Cell complexes. See Hatcher [T5i for definitions and basic facts about cell complexes. We 
write Xi for the set of all j-dimensional cells of X, and Xu\ instead of Hatcher's X 1 for the i- 
skclcton Xq U Xi U • • • U Xi. The notation X d indicates that X has dimension d. We borrow some 
standard terminology from the theory of simplicial complexes: a cell of X not contained in the 
closure of any other cell is called a facet of X, and we say that X is pure if all its facets have the 
same dimension. 

Let R be a ring (if unspecified, assumed to be Z), and let Ci(X) denote the i th cellular chain 
group of X, i.e., the free i?-module with basis {[F]: F G Xi}. We then have cellular boundary and 
coboundary maps 

d x y. d(X) -> Q-^X), 
d* x y. d-^X) ^ d(X), 

where we have identified cochains with chains via the natural inner product. We will abbreviate the 
subscripts in the notation for boundaries and coboundaries whenever no ambiguity can arise. 

We will often regard d i (resp. d*) as a matrix whose columns and rows (resp. rows and columns) 
are indexed by Xi and Xi_\ respectively. 

The i th (reduced) homology group of X is Hi(X) — ker(9 i )/ im(d i+1 ), and the i th (reduced) Betti 
number fti(X) is the rank of the largest free i?-module summand of Hi(X). 

2.2. Laplacians. For < i < dimX, define linear operators L x i i , L x A i , L x \ on the vector space 
Q(X) by 

Lx,i = ^i+i^f+i (the up-down Laplacian), 

^Xi = (the down-up Laplacian), 

l}xi = + Lj? ti (the total Laplacian). 

Let s" d (X), sf u (X) and s' ot (X) denote the spectra of the respective Laplacians, that is, the 
multisets of their eigenvalues. Each of these multisets has cardinality \Xi\ (counting multiplicities), 
because the Laplacians can be represented by symmetric square matrices of that size. Instead of 
working with multisets, it is often convenient to record the eigenvalues by the generating functions 

E'(X;r)= J2 ^ 
\es'(X) 

dim X 

E-(X;r,t)= ]T J2 tV " 

where • € {ud, du, tot}. 

The various Laplacian spectra are related by the identities 

Ef(X;r) = Ef^(X;r), (3a) 

E?*(X;r) ^ Ef(X;r) + Ef»(X;r) = Ef(X;r)+Ef\(X;r), (3b) 

Ef{X; r) = E**{X; r) - E^(X; r), (3c) 

eqn. (3.6)], where / = g means that the two polynomials /, g are equal up to their constant 
coefficients. It follows that 

i 

Ef^(X-r) ± Ef(X-r) ± £(-1)-^° \X; r). (4) 
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Therefore, each of the three families of generating funtions 

{E\°\X-r): 0<i< dimX}, {Ef d (X;r): 0<i< dimX}, {Ef u {X;r): 0<i< dimX} 

determines the other two. (The constant coefficients, which represent the multiplicity of the zero 
eigenvalue, can always be found by observing that Ef d (X, 1) = E? U (X, 1) = E\ ot {X, l) = \X t \.) 

The bookkeeping differs slightly from simplicial complexes: each cell in a cell complex has di- 
mension > 0, whereas every simplicial complex includes the (— l)-dimensional face 0. Therefore, for 
instance, E du (X,q) = for all X. 

2.3. Product and sum formulas for the total Laplacian. For this section, it is convenient to 
regard the boundary map of a cell complex X as a linear endomorphism d x : C(X) — > C(X), where 
C(X) — i Ci(X). That is, d x = J2 i d x t ; equivalently, for each i, d x • is the restriction of d x to 
the group Ci(X) of cellular i-chains. 

If X and Y are disjoint cell complexes, then 

E tot (X U Y; r, t) = E tot (X; r, t) + E tot (Y; r, t), 

because the k th total boundary map (resp., Laplacian) of X U Y is just the direct sum of the 
corresponding boundary maps (resp., Laplacians) of X and Y. 
Let Z = X x Y. Then the set of fc-cells of Z is 

k 

Zk = [J Xi x Yfc_i, 

i=0 

and on the level of chain groups we have 

k 

C k (Z) = ®Ci(X)®C k -i(Y). 

The boundary map d z is defined as follows. Let / 6 Xi and g £ Then (/, g) e Z k , and 

the boundary map acts on the corresponding cellular chain [/] <g> [g] e Ck(Z) by dim/ = i and 
dim g = k — i, 

dzdf] ® [g]) = d x [f] ® [g] + (-l) dim/ [/] «g) d Y [g\. 
More simply, we may write 

k k 

d z = ^]9 x (g)id+(-l)Md®5 y , d z = ^d x id®d^. 

j=0 i=0 

From this we calculate 

Lf = d z d z = (-iy- 1 d x ®d^+L u x i ®id + id®L¥ i + (-iyd x ®d Y , 

Lf = d z d z = (-iyd x ®d Y + Lf®id+id®Lp + (-iy- 1 d x ®d Y , 

= L % x ® id + id ®L^\ 
That is, the matrix for is block-diagonal, with blocks 

L%®id + id®I%* k , L x \ ® id + id ®L^ k _ 1 , L^ k ® id + id^L^. 

Therefore, 

4 ot (Z) = {A + M : A e s * ot (x), M G 4 ot (r), i + j = fc} 

as multisets. From this we obtain the following product formula. 
Theorem 2.1. Let X,Y be cell complexes. Then 

E tot (X x Y; r, t) = E tot (X; r, t)E tot (Y; r, t). 
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2.4. Spanning trees. 

Definition 2.2. Let X d be a cell complex, and let k < d. A cellular spanning k-tree (for short, 
CST or fc-CST) of A is a fc-dimensional subcomplex Y" C A such that Y^-i) = ^(k— l) an d 

# fc (Y)=0, (5a) 

l-fffc-iPOl < co, and (5b) 

in! = 1^1 -yfffcC^) + ^ fe -i(J ! «:). (5C) 

We write Tk(X) for the set of all fc-CST's of X, sometimes omitting the subscript if k = d. Note 
that T k (X) = Tk{X {j) ) for all j > k. 

A O-dimensional CST is just a vertex of A. If A is a connected 1-dimensional cell complex — that 
is, a connected graph — then the definition of 1-CST coincides with the usual definition of spanning 
tree: a spanning subgraph of A that is acyclic, connected, and has one fewer edge than the number 
of vertices in A. 

Proposition 2.3. Let Y C A be a k- dimensional subcomplex with Yffe-i) = X( k _iy Then any two 
of the conditions (|5aj) . (|5b|) . ((5c|) together imply the third. 

The proof is identical to that of [8j Prop. 3.5]. 

Definition 2.4. A cell complex A is acyclic in positive codimension, or APC for short, if (3j(X) = 
for all j < dim A. 

As in the setting of simplicial complexes, APC-ness is the "right" analogue of connectedness for 
graphs, in the following sense. 

Proposition 2.5. A cell complex has a cellular spanning tree if and only if it is APC. 
The proof is identical to that of [8j Prop. 3.7]. 

2.5. The Cellular Matrix- Tree Theorem. Throughout this section, let X d be an APC cell 
complex with d > 1. Define 

T k =r k (X)= Yl \Hk-i(Y)\ 2 forO<fc<d, 

Y£T k (X) 

7T fc = 7r fe (A) = J"! A for 1 < fc < 

Observe that r = |A |, because a O-dimensional CST is just a vertex. In addition, wc define 

7r =7ro(X) = |X |, 7T_! = 7T_ 1 (A) = 1. 

While it might seem more consistent to define ttq — 1 (because sH^A) = 0), several subsequent 
enumeration formulas (such as Corollarv l3.5l and Theorem l6.2[) can be stated much more conveniently 
with this convention for 7ro . (By way of motivation, in the case of a simplicial complex on n vertices 
(see [8]), it is in fact true that ttq = n, because of the presence of a (— l)-dimcnsional face; the 
Laplacian is is not a void matrix, but rather the lxl matrix with single entry n.) 

Abbreviate^ = (3i{X) and d — d Xd . Let T be a set of d-cells of A of cardinality \X d \— Pd+Pd-i — 
I Ad | — (3d, and let S be a set of (d — l)-cells such that |5| = \T\. Define 

X T = T U A(d-i), S = A( rf _!) \ S, Xg — S U X( d _ 2 ), 
and let d s T be the square submatrix of d with rows indexed by S and columns indexed by T. 

Proposition 2.6. The matrix d ST is nonsingular if and only if Xt £ 7d(X) and Xg G 7d-i(A). 
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Proof. We may regard d s T as the top boundary map of the ci-dimensional relative complex T = 
(Xt,X s ). So d ST is nonsingular if and only if Hd(Y) = 0. Consider the long exact sequence 

-»• H d {X s ) -> H d (X T ) -> H d (Y) -> H d -i(X § ) -> # d _i(X T ) -> H d _i(r) (6) 

If Hd(T) ^ 0, then Hd(Xx) and Hd-i(X s ) cannot both be zero. This proves the "only if" direction. 

If ff rf (r) = 0, then # d (A 5 ) = ( since dim Ag = d-1), so © implies H d {X T ) = 0. Therefore A T 
is a d-tree, because it has the correct number of facets. Hence Hd-ii^r) is finite. Then ([BJ implies 
that Hd-i(Xg) is finite. In fact, it is zero because the top homology group of any complex must be 
torsion-free. Meanwhile, X s has the correct number of facets for a (d — 1)-CST of X, proving the 
"if" direction. □ 

Proposition 2.7. Ifd ST is nonsingular, then 

ldetd , _ \H d -i(X T )\ ■ \H d - 2 (Xs)\ _ \H d -i(X T )\ • \H d -2(X s )\ 
S ' \H d -2(X T )\ Wd-2{X)\ 

Proof. As before, we interpret d s T as the boundary map of the relative complex Y — (Xt, Xg). So 
d s T is a map from Zl T l to Zl T l, and l) T \ / d s T {l) T \) is a finite abelian group of order | det d ST \. On 
the other hand, since Y has no faces of dimension < d — 2, its lower boundary maps are all zero, so 
| det d s T \ = \Hd-i(Y)\. Since Hd~2(Xr) is finite, the desired result now follows from the piece 

-> H d -i(X T ) -> ^_ x (r) ->■ ff d _ 2 (A 5 ) -> ff d _ 2 (A T ) -> (7) 
of the long exact sequence (O ■ □ 

We now can state our main result connecting cellular spanning tree enumeration with Laplacian 
eigenvalues. 

Theorem 2.8 (Cellular Matrix- Tree Theorem). Let d > 1 and let X d be an APC cell complex. 
Then: 

(1) We have 

TdiX^d-^X) 

ttAX) = . 

(2) Suppose that d > 0. Let L = L^j 1( let U be the set of facets of a (d — \)-CST of X, 
and let Ljj denote the reduced Laplacian obtained by deleting the rows and columns of L 
corresponding to U. Then 

r d (X)= fo-sWI 2 det L V . 
dK ' \H d -2(Xu)\ 2 

Proof of Theorem \2.8[ (1). The Laplacian L is a by |X<j_i| square matrix. Since X is APC, 

we have rankL = \X d \ — (3d — \Xd\ — $d + Pd-i- Let x{L;y) — det(yl — L) be the characteristic 
polynomial of L (where / is an identity matrix of the same size), so that TTd(X) is given up to sign 
by the coefficient of y\ x *-i\-\ x <i\+P<' i n x{L;y). Equivalently, 

nd = ^3 detLu = det Lu (8) 

5CX d _i SCXd-i 
|S|=rankL \ S \ = \X d \-fJ d 

where U = Xd-i\S in each summand. By the Binet-Cauchy formula, we have 

det Lu = ^2(detd s ^ T )(detd* S T ) = ^ (det<9 S T ) 2 . (9) 

TCX d TCX d 
\T\ = \S\ \T\ = \S\ 
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Combining ([8} and (j9]), applying Proposition 12. 6[ and interchanging the sums, we obtain 

n d = E ( detd s.T) 2 

T:X T &T d (X) S:X s (LT d -i(X) 

and now applying Proposition 12 . 71 yields 



— E E 

T:X T £T d (X) S:X s eT d -i(X) 



\H d ^(X T )\-\H d ^ 2 (Xs)\ \ 2 
\H d - 2 (X)\ J 



£ \H d ^(X T )\^ Yl \H d -2{X s )\ 

K T:X T dT d (X) J \S:X s eT d - 1 (X) 

~ \Hd-2(X)\ 2 
as desired. □ 



In order to prove the "reduced Laplacian" part of Theorem l2.8[ we first need the following lemma. 

Lemma 2.9. Let U be the set of facets of a (d - l)-CST of X, and let S = X d -i\U. Then 
\S\ = \X d \ — f3 d (X), the number of facets of a d-CST of X . 

Proof. Let Y = X/ d _iy, in particular, 

\Y t \ = \X t \ for ^< d-1 and (3 e (Y) = f3 e (X) for I < d - 2. (10) 

Therefore, \U\ = |Y d _i| - fa-^Y) + /3 d _ 2 (F) = \X d ^\ - /3 d _i(Y), so \S\ = fl,_i(Y) by Proposi- 
tion 12.31 Meanwhile, the Euler characteristics of X and Y are 

d d 

x (x) = = £(-i)*ft(x), 

i=0 i=0 
d-1 d-1 

x(^) = E(- 1 ) i i^i = E(- 1 ) i ft( F )- 

i=0 i=0 

By (fTOl) . we see that 



= (-1)1^1 = (-l^W + t-lJ^^lW-t-l)^^-!^) 

from which we obtain \X d \ = (3 d (X) - fa-x{X) + /3d-i(Y). Since X is APC, we have /3 d -i(X) = 0, 
so \S\ = p d -i(Y) = \X d \ - (3 d {X) as desired. □ 

Proof of Theorem ] 2. 8\ (2). By the Binet-Cauchy formula, we have 

detL[/= ( detd s,T)( detd *s,T) = E ( detd s,T) 2 - 

T: \T\ = \S\ T: |T| = |S| 

By Lemma 12.91 and Proposition 12.61 d ST is nonsingular exactly when Xt £ T d (X). Hence 
Proposition 12.71 gives 

f \H d - X {X T )\ ■ \H d _ 2 (Xu)\\ 2 



det Lu = ^ 

T:X T er ti (x) 



l^-2(X)| 



= 7^ T^TTT" l^d-i(A T )| = -y~: — — r d (AJ, 

l#d- 2 (*)l 2 T:XTeTd(x) I#d-2p0l 2 

which is equivalent to the desired formula. □ 
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We will often work with complexes that are in fact Z-acyclic in positive codimension (as opposed 
to merely Q- acyclic ) , and whose Laplacians have nice forms. In this case, the following formula is 
often the most convenient way to obtain tree enumerators from Laplacian eigenvalues. 

Corollary 2.10. Let X d be a cell complex such that Hi(X,Z) = for all i < d. Then, for every 
k < d, we have 

k 

n (_!)*-< 

i=0 

Proof. Applying Theorem 12.81 fl) repeatedly, we obtain 

/ k \ k 

TTfc 7Tfc / -p-r (_!)*-* \ _ T-r (~l) fc - 



Tk = = Tk-2 

Tfc-1 TTk-1 



\i=l / i=0 

For later use (when we study duality in Section [5]) , define 

Lj k = uJk{X) = Yl A. 

o^Aes* ot (Jf) 

Note that equation (|3bp implies that uj^ = "Kh^k+i for k > 0. Moreover, luq = tt\. Solving for the 
7r's in terms of the w's gives 

j-i 

n(_i)*-j+i 
u j 

3=0 

and substituting this formula into Corollary 12.101 gives 



Tfe = 11 IK = 11 11 % = 11^ 



i=0]=0 ]=0i=] + l ]=0 

l[^ h - i+1 ^. (11) 
3=0 



2.6. The Weighted Cellular Matrix- Tree Theorem. As before, let d > 1 and let X d be a cell 
complex that is APC. Introduce an indeterminate Xf for each / G X dl and let Xf = x 2 ?. For every 

T C A^, let xt = 11/iet x f ano - ^ = x t- To construct the weighted boundary matrix dx,d, we 
multiply each column of d x d by Xf, where / is the corresponding d-cell of X . We can accordingly 
define weighted versions of the coboundary maps, Laplacians, etc., of Section |2~51 as well as of the 
invariants and r^. We will notate each weighted object by placing a hat over the symbol for the 
corresponding unweighted quantity. Thus 7Tfc is the product of the nonzero eigenvalues of L x i k _ 1 
(for k > 1), and 

f fe =f fe (A)= \H k -i(Y)\ 2 X Y . 

YET k (X) 

Meanwhile, any unweighted quantity can be recovered from its weighted analogue simply by setting 
Xf = 1 for all / 6 Xd. 

Proposition 2.11. Let T C X d and S C A rf _i, with \T\ = \S\ = \X d \ - /3 d . Then detd s , T = 
xt det d s T is nonzero if and only if Xt G T d {X) and Xg G 7d-i(A). In that case, 

±det§ = \H d ^X T )\-\H d - 2 (X s )\ = \H d -i(X T )\ ■ \Hd-2{X§)\ 
\H d . 2 (X T )\ \H d -2(X)\ 
Proof. The first claim follows from Proposition ^. 6[ and the second from Proposition 12. 71 □ 

It is now straightforward to adapt the proofs of both parts of Theorem l2.8l to the weighted setting. 
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Theorem 2.12. Let d>l, let X d be an APC cell complex, and let L = L x i d _ 1 . 

(1) We have 

= UX)r d ^X) 
dK 1 \H d - 2 (X)\ 2 

(2) Let U be the set of facets of a (d— 1)-CST of X , and let L\j be the reduced Laplacian obtained 
by deleting the rows and columns of L corresponding to U. Then 

f d (X)= fo- 2W|2 detLu. 
dK ' \H d - 2 {Xu)\* 

3. Cubical complexes 

We now specialize from arbitrary cell complexes to cubes and cubical complexes. We retain the 
notation and terminology of Section [2] for cell complexes. 

The n-cube Q n is the face poset of the geometric n-cube Q n , the convex hull of the 2" points 
in 1" whose coordinates are all or 1. We will usually blur the distinction between the polytope 
Q n and its face poset Q n . The faces of Q n are ordered n-tuples / = (fi, ■ ■ ■ , f n ), where /j G 
{0, 1, ©}. Intuitively, we think of © as the open interval (0, 1). For example, the cell (0, ©, 1, ©) G Q4 
corresponds to the 2-cell 

{0} x (0,1) x {1} x (0,1) C Q A C M 4 . 

The order relation in Q n is as follows: / < g iff /, < gi for all i G [n], where < ©; 1 < ©; and 0, 1 
are incomparable. (If necessary, we can regard Q n as containing a unique minimal element 0, with 
undefined direction and dimension —1.) 

The direction of a face / is defined as dir(/) = {i G [n] : fi = ©}, so that dim(/) = |dir(/)|. 
Notice that dir(/) C dir(g) whenever f < g, although the converse is not true. The poset Q n is 
ranked, with (™)2™ -i faces of rank i for < i < n. 

A proper cubical complex X is an order ideal in Q n . This is a combinatorial object with a natural 
geometric realization X as the union of the corresponding faces of the polytope Q n . Note that this 
is a much more restrictive definition of "cubical complex" than as simply a cell complex all of whose 
faces are combinatorial cubes. The reason for working with this smaller class of cubical complexes 
is that the cells and boundary maps can be described combinatorially, as we now explain. 

Let / and g be faces of Q n of dimensions i — 1 and i respectively. If / < g, then we may write 
dir((/) = {01, . . . , cii}, dir(/) = {01, . . . , a}, . . . , ai}, with both direction sets listed in increasing order. 
Then the relative orientation of the pair /, g is 

Vl) J if /«.,=(), 
if/,, = l. 

Meanwhile, if / j£ g, then we set e(f, g) — 0. 

Let R be a coefficient ring (typically Z or a field), and let Ci(X) = Ci(X,R) be the free abelian 
group on generators [/] for / G X j. The i th cubical boundary map of X is 

3 x y. Ci(X) ^ Ci-i(X) 

and the i th cubical coboundary map of X is 

d* x y. Ci-!{X) ^ d(X) 



£(f,g) 
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As before, we define the i up-down, down-up, and total cubical Laplacians as respectively 



t ud 



- du 



n* r du a* o t tot run . r a 

J X,i+l°X,i+l> ^X,i — a X,i°X,ii ^X,i — ^X,i + X, 

The map d x i is in fact the cellular boundary map of Q n as a cell complex; see [TTJ §4] . So the 
techniques of Section [5] can be applied to count cubical spanning trees. 

Example 3.1. A fundamental example is the complex A = Q 1: with A — {0, 1}, X\ = {©}. The 
boundary, coboundary and Laplacian matrices of X are 

T 



1 



j ud 
^0 



T tot 

L Q - 



1 



-1 



-1 1 



t l r x 



L tot 



1 



tr l 



(13) 



^ = 

0?=[-l 1], 
and so its spectrum polynomial is 

dimX 

E^(Q i; r,t) = Yl E 

»=o Aes* ot (x) 

3.1. Prisms. We now consider the important prism operation, which is the cubical analogue of 
coning a simplicial complex. 

Definition 3.2. Let X C Q n be a proper cubical complex. The prism over X is the subcomplex of 
Qn+i defined by 

= {/ = (A, • - • , U+i) e Q„+i : (A, ..,/„)£!}. 

Note that 

(PX)i = {(fx, ...,/„, 0), (A, ...,/„, 1) : (/i, ...,/„)£ AJ 
U{/ = (/i,...,/ nj ®): (fx,...,fn) 

As a cell complex, PA is just the product X X Qj. Therefore, we can use the product formula, 
Theorem 12. 11 to write down the Laplacian spectra of PA in terms of those of A. The formula can 
be stated in several equivalent ways, all of which will be useful in different contexts. 

Theorem 3.3. Let X d be a proper cubical complex. Then 



£ tot (Pl; r, t) = (1 + r 2 + tr 2 )E tot (X; r, t), (14a) 

El ot (PX; r) = (1 + r 2 )El ot (X; r) + r 2 E\^ x (X; r), (14b) 

s* ot (PA) = {A: Aesf(I)} U {A + 2 : A £ s* ot (A)} U {^ + 2:^^(1)}, (14c) 

sf(PX)^{\:\esf(X)} U {A + 2: Ae^W}, (14d) 

for all < i < d, where U denotes multiset union. 

In particular, the prism operation preserves Laplacian integrality. 

Proof. Equation (|14al) follows from Theorem 12.11 together with equation (|13p . and (|14bl) and (|14c|) 
are just rephrasings of (|14a[) . 

To prove (| 14d[) . we proceed by induction on i. For i = 0, the formula follows from (j!4c|) . because 
sJj d (A) = 4 ot (A) and s^KX) = 0. For i > 1, we have 

sf(PX)^s? t (PX)\s^ 1 (PX) 

= s™(PX) \ (s^(X) U {A + 2: A e S «(I))) 

= s 4 ud W U {X + 2: Xesl ot (X)} U s^(A) U {A + 2 : A e (A)} 

\( S ^(A) U {X + 2:\est\(X)}) 
= sf(X) U {A + 2: Ae4°'(I)}. 

□ 
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3.2. Laplacian spectra of cubes. As a consequence of Theorem I3.3[ we obtain a formula for the 
Laplacian eigenvalues of Q n , and thus for the torsion-weighted spanning tree enumerators Tk{Q n )- 

Theorem 3.4. Cubes and their skeletons are Laplacian integral. Specifically, for all n > 1, we have 



E tot (Q n ,r,t) = (l + r 2 +tr 2 r = £ ** r 2k (1 + r 2 ) 



fe=0 



„2\n-k 



n-l 



E ud (Q n ;r,t) = £ t* 



k=0 



E 

3 =k+i 



n\ ( j - 1 



.2,/ 



(15) 



(16) 



Proof. The formula for Ej. ot (Q n ; r, t) follows from Theorem 13.31 since Q n can be identified with the 
n-fold product Qi X • • • X Qi as cubical complexes (indeed, as cell complexes). 

By |5c]). we can obtain E ud (Q n ; r, t) from E tot (Q n ; r, t) by deleting all the r-free terms (i.e., those 
which correspond to zero eigenvalues) and dividing by 1 + 1. The only such term is a single 1 (from 
the k = summand). Therefore, 

(l + r 2 (l + t))"-l 



E ua (Q n ;r,t) = 



l + t 



r 2j (l + t) j 



l + t 



E 



i+t 

^(i+tr 1 = E( - PE 



fe=0 



i-i 



E* fe 

fc=0 



E 

j=fe+i 



n \ a - 1 



,2.; 



□ 



Corollary 3.5. Let n > 1 and 1 < fc < n. TTien 



r fc (Q n ) = [] (2j)(")^ : 
j=fc+i 



Proof. Theorem 13.41 implies that 



^(Q«)-n^) a)(u) =n^) a)a=;) 

for < i < n (adopting the convention that (?) = for a < b). Applying the alternating product 
formula, Corollary 12.101 gives 

{ _ 1)k -i 



Tk 



(Qn) = 



(-l) k - 



2 „(-l) fc -Q 



n(2i)(-)«-) 



= 2"(- 1 ) fe IJ(2j)(")(^(- 1 ) fc "(ti)). 

The j — 1 factor in this product simplifies to 2"'~ 1 ) , canceling the initial factor of 2™(~ 1 - ) . For 
2 < j < fc, the sum in the exponent vanishes, and for k + 1 < j < n, it simplifies to (i^i) (this can 
be seen by repeatedly applying Pascal's identity), giving the desired formula. □ 
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4. Weighted Laplacians of cubes 



In this section, we study a weighting that associates a Laurent monomial to each face of Q n , 
giving finer information about Laplacian spectra of cubes. 

4.1. Algebraically weighted eigenvalues. Let X C Q n be a proper cubical complex, and intro- 
duce commuting indeterminates Xi, yi, qi for i 6 N. Weight each face / = (/i, . . . , /„) G Q n by the 
monomial 

& = n qi n x n y *- 

i:fi=® i:/i=Q i:/«=l 

Define the algebraically weighted cubical boundary map by 

dx,k- C k {X) ^C k -i(X) 



E 



(17) 



so that the corresponding weighted coboundary map is 

§* x y. CVi(X) C*(X) 



seA fc ?/ 



(18) 



It is easy to check that dd = d*d* — 0. (This vital equality would fail if we had defined the weighted 
boundary more "combinatorially naturally" by [g] t-¥ X)/eA fc x £ (/i flO£s[/]0 

The /c" 1 up-down, down-up, and tota^ algebraically weighted cubical Laplacians are respectively 



fc+1; 



r tot rud 

^X,fe — ^X,k 



J X,k- 



As in the unweighted case, it is convenient to record the eigenvalues as generating functions: 



E' k (X;r) = 



dim AT 



E ^\ 



E'(X;r,t)= £ 

k=o Aesj(x) 

where • £ {ud, du, tot}. 

Example 4.1. Consider the complex X — Q\ (see Example 13. 1| . whose edge we regard as lying in 
direction i. The vertices have weights Xi and j/j, and the edge has weight qi. The weighted boundary, 
coboundary and Laplacian matrices are thus 



ft 



-qi/x. 
_ Qi/Vi 

91 = [-Q.il x% Qi/Vi] , 
and so the weighted spectrum polynomial is 

E tot (Qv,r,t) = Y, 

i 

where 



f ud 

Li 



o 

du 



E 

Aes* ot (Qj 



7- tot 
^0 



if' = [ 



t l r A = 1 



Qi/x 2 i -qf/xzVi 

-qf/xiVi qf/yf 
qf/xf + qf/yf] , 

+ tr u ' 



(19) 
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Just as in the unweighted setting, the total eigenvalue generating function E tot (X;r, t) is multi- 
plicative on products of cubical complexes, that is, 

E tot {X x Y-r,t) = E tot (X;r,t)E tot (Y;r 7 t). (20) 

This formula is proved in exactly the same way as its unweighted analogue Theorem 12.11 (which can 
be recovered from (|2"0"| by setting % = Xi = = 1). In particular, if X C Q n -\ is a proper cubical 
complex, then the prism PX C Q n is just the product of X with a copy of Qi lying in direction n. 
Hence 

E tot (PX;r,t) = (1 + r u " + tr u ") E tot (X;r,t). (21) 
Note that specializing qi = = y.i = 1 makes Uj = 2, and so recovers the first assertion of 
Theorem 13.31 The corresponding recurrence for the multisets of eigenvalues is 

sT'iPX) = {X: \e§l ot (X)} U {X + u n : Ae^(I)} U { f i + u n : f i e st\(X)} . (22) 

We now calculate the eigenvalues of the weighted Laplacians of the full cube Q n - 

Theorem 4.2. For all n > 1, we have 



E tot (Q n ;r,t) = Y[(l + r u * +tr Uk ) = £V 



fc=i 



i=0 Lfe=i 



E 



f J A: 



^ ud (Q„;r,t) = ^(i + tr 1 e J = X) 4< 

j=l i=0 



E 

j=i+l 



(23) 
(24) 



where ej denotes the j elementary symmetric function in the forms r Ul , . . . , r Un . 

Proof. The first equality of (f2"3")l follows from iterating the weighted product formula (f2"0"|) . Extracting 
the coefficient to find E t i ot (Q n - 1 r) and abbreviating ua — SigA u i f° r & C [n] we have 



E^ E 



AC[n] ye A 

Ul=« 



EE E< 

fc=i CC[n] ACC 
\C\=k \A\=i 



AC[n] 
\A\=i 



BC[n]\A 



E E 

k=i CC[n] 
C|=fc 



E 

k—i 



which is the second equality of (|2"3")l . 

Finally, by the weighted analogues of the equalities (|3c|. we can obtain E ud (Q n ;r,t) from 
E tot (Q n ',r,t) by deleting all the r-free terms (i.e., those which correspond to zero eigenvalues) and 
dividing by 1 + t: 

t f[ 1 + (1 + t)r" fc j - 1 



E^{Q n ;r,t) 



1+t 

n n j—1 

j2n + ty-^ = j:j: 

j=l J=l i=0 



i-i 

i 



t*e,- 



E^ 

i=0 



E 



i-i 
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□ 



Extracting the t % coefficient from formula ([24)) of Theorem 14.21 gives 



Er(Qn,r) 



AC[n] 
\A\>i+l 



That is, the nonzero eigenvalues of L" d (Q n ) are the expressions ua, each occurring with multiplicity 
('^'r 1 ). Therefore, the product of nonzero eigenvalues — that is, the algebraically weighted analogue 
of the invariant 7Tj(Q n ) — is 



n 



Aesr(Q„) 

A/0 



n 

AC[n] 
\A\>i+l 



n 

AC[n] 

\A\>i+l 




(25) 



4.2. Weighted tree enumeration. We now consider the polynomial generating function 



E 

YeT k (X) g£Y k 



for a proper cubical complex X. As usual, the main case of interest is X = Q n . In principle, the 
invariants Tk{Q n ) can be computed in terms of the eigenvalues of weighted Laplacians, using The- 
orem [2TT2] Those "combinatorially weighted" Laplacians look similar to (in fact, simpler than) the 
algebraically weighted Laplacians discussed in the previous section, but they do not come from well- 
defined boundary maps (i.e., whose square is zero) and their eigenvalues are not even polynomials. 
On the other hand, there is strong evidence for the following formula. 



Conjecture 4.3. 



Tk(Qn) = (?1 ■■■Qn) 



n-1 

E 

! = fc-l 



n 

AC[n] 
\A\>k+l 




jeA\i 



\A\-2 

k - 1 



The conjectured formula is similar to equation (|25[) : specifically, clearing denominators from the 
square-bracketed expression in (|25[) indexed by A C [n] gives the corresponding factor indexed by A 
in the right-hand side of the conjecture. Our original goal in proving formulas such as Theorem 14.21 
was to prove Conjecture 14.31 by translating between the algebraically and combinatorially weighted 
Laplacians; this approach had succeeded in the case of shifted simplicial complexes [5], but it is not 
clear how to do that here. 

Conjecture 14.31 can be verified computationally for small values of n and k. The case k = 1 is 
Theorem 3 of [20], and the case k = n — 1 can be checked directly, because an (n — l)-spanning tree 
of Q n is just a subcomplex generated by all but one of its (n — l)-faces. 

We suggest a possible avenue for proving Conjecture 14.31 First, note that the weighted spanning 
tree enumerator is homogeneous in the variable sets {xi,yi,qi},...,{x n ,y n ,q n }, so we lose no 
information by setting q± = ■ ■ ■ = q n = 1 on the right-hand side, obtaining the simpler formula 

'\A\-2 
. k-1 



F(n,k)= J] 



AC[n] 
A|>fc+1 



E 

ieA 



jGA\i 



Second, we observe that F(n, k) satisfies the recurrence 

F[n,k) = G([n],n- l,k)F[n,n- 



CELLULAR SPANNING TREES AND LAPLACIANS OF CUBICAL COMPLEXES 



17 



where 

G(S,a,b)= J] F(A,b). 

ACS,\A\=a 

5. Shifted cubical complexes 

Inspired by the notion of shifted simplicial complexes (see, e.g., [17j [7]), we define the class 
of shifted cubical complexes. These shifted cubical complexes share many properties with their 
simplicial counterparts. In particular, they are Laplacian integral, constructible from a few basic 
operations, and arise as order ideals with respect to a natural relation on direction sets. The guiding 
principle is that directions in cubical complexes are analogous to vertices in simplicial complexes. 

We first generalize our definitions slightly to be able to work with cubical complexes on arbitrary 
direction sets. That is, if D is any set of positive integers, a cubical complex with direction set D is a 
family A of ordered tuples (fi)ieD, where /, G {0, 1, ©} for all i, closed under replacing ®'s with 0's 
or with l's. The direction and dimension of faces are defined as before: dir(/) = {i G D : fi = ©}, 
and dim / = | dir(/)|. We will frequently need to regard the direction of a face as a list in increasing 
order; in this case we write dir(/)< instead of dir(/). 

Definition 5.1. A shifted cubical complex is a proper cubical complex A C Q n that satisfies the 
following conditions for every /, g G {0, 1, ©}™ with dim / = dim g: 

(1) X contains the full 1-skeleton of Q n . 

(2) If / G X and dir(/) = dir(g), then g G X. 

(3) If g G X and dir(/) < precedes dir(g) < in component- wise partial order (that is, the i th 
smallest element of dir(/)< is less than or equal to that of dir(g) < for every i), then / G X. 

The first condition is analogous to the requirement that a simplicial complex "on vertex set V" 
actually contain each member of V as a vertex. In light of condition ([2]), it would be equivalent 
to replace ([T]) with the condition that X contains at least one edge in every possible direction. 
The second condition reflects a symmetry between the digits and 1, while the last condition is the 
cubical analogue of the definition of a shifted simplicial complex (as a complex whose faces, regarded 
as collections of vertices, form an order ideal in component- wise partial order). 

5.1. Near-prisms. Bjorner and Kalai [3] introduced the concept of near-cones, computed their 
homotopy type, and showed that shifted simplicial complexes are near-cones. The cubical counter- 
part of a near-cone is a near-prism. We develop the basic facts about near-prisms in parallel to the 
presentation of Section 5] , in order to prove that shifted cubical complexes are Laplacian integral. 

Throughout this section, let X C {0,1,©}™ be a d-dimensional proper cubical complex. For a 
direction i G [n], define the deletion and link with respect to i as follows: 

del i X = {/\/ i : /el, ft^®}, 

link i X = {f\f i : /eX, /»• = ©}. 

The deletion and link in direction i are proper cubical complexes with direction set [n]\i. Mean- 
while, given a complex X on direction set [n]\i, define the prism of X in direction i as follows: 

P,A = {(/x, ...,/„) G {0, 1, ©}" : f\fi G X}. 

If X C Qn-i, then the prism PA defined in Section l3~T1 is P„A in this notation. 

Observe that P;A naturally contains two isomorphic copies of X: one consisting of all faces of X 
with inserted in the i th digit, and one consisting of all faces of X with 1 inserted in the i th digit. 
We denote these subcomplexes by l X and VX respectively. 

Definition 5.2. A cubical complex A is a near-prism in direction i if (1) the boundary of every 
face of deb A is contained in link^ A, and (2) both □'(deb A) and l*(deb A) are contained in A. 
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If X is a near-prism in direction i, then it admits the decomposition 

X = Q^delj X) U l^delj X) U P^link, X). (26) 
The following fact is not difficult to prove directly from the definition of a shifted cubical complex. 

Lemma 5.3. X is a shifted cubical complex iff X is a near-prism in direction 1 and both linki X 
and deli X are shifted cubical complexes with respect to the direction set {2, . . . , n}. 

In particular, shifted cubical complexes are iterated near-prisms. As in [7], this characterization 
will help describe the Laplacian eigenvalues of a shifted cubical complex. 

First, we introduce notation to work with weakly decreasing sets of nonnegative integers, defined 
up to =-equivalence (i.e., with an indeterminate number of trailing zeroes). The symbol 2™ will 
denote the sequence (2,2,..., 2) of length m. If a = (ai > ct2 > . . .) and b = (61 > 62 > ■ ■ •) are 
two sequences, then a + b denotes the sequence (ai + 61 > 02 + 62 > • • ■)• I n particular, 2™ + a is 
the sequence derived from a by adding 2 to each of the first m entries of a, padding the end of a 
with O's if necessary. For instance, 2 8 + (7, 5, 5, 2, 1) = (9, 7, 7, 4, 3, 2, 2, 2). 

Let s(X) denote the sequence of non-zero eigenvalues of L^ d _ 1 , written in weakly decreasing 
order. 

Theorem 5.4. Let X d be a pure cubical complex which is a near-prism in direction 1, and let £ 
denote the number of facets 0/ linki X. If X is also a prism, then 

s(X) = 2 £ + s(linkiX). 

Otherwise, 

s(X) = s(deli X) U (2 £ + (s(deliX)Us(linkiX))). 

Proof. First, suppose that X is not a prism. Then dim(deli X) = d and dim(linki X) = d — 1. Let 
X' = deli X. Since X is a pure near-prism, we see that 

X' d = deb X and X' d _ t = linki X. 

So (PiX') d = 1 deli X U l 1 deb X U Pi linki X = X by equation ^ above. Therefore 

sf^X) = sf^PX') 

= sti(x') u (2W- 1 i+4°! 1 (x')) 

= s d d _i(X') U (2\ x '^\+(sti(X') U a£_ 2 (X'))) (27) 

= S ^i(deliX) U (2K linklX )- 1 l + ( S ^ 1 (del 1 X)U S ^ 2 (linkiX))), 

where the second line uses (|14d|) . and all integer sequences are listed in weakly decreasing order. 

The proof is similar if X is a prism. In this case, X = PX' where dimX' = d — 1, and 
X 1 = linki X, so again X' d _ 1 = linki X. Equation (|27p then applies again, but s d '^_ 1 (X') consists of 
all O's, because dimX' = d— 1. Therefore 

sfn{X) ^l^ 1 ^- 1 ' +s3l 2 (linki X). 

□ 

Corollary 5.5. Shifted cubical complexes are Laplacian integral. 

Proof. We proceed by induction on the number of directions. If a shifted cubical complex X has 
only one direction, then X = Qi, which is Laplacian integral by Theorem 13.41 

Now let X be a ci-dimensional shifted cubical complex on more than one direction. Let Y be the 
j- dimensional pure skeleton of X, i.e., the subcomplex consisting of the j-dimensional faces of X and 
all of their subfaces. It is clear that s^^X) = ^ U (X) = s^ u (y) = sj£ x (y), since s^ u depends only 
on Xj and Xj-i. Therefore, we only need to show that Y is Laplacian integral, but Y is pure by 
definition, so we may apply Theorem 15.41 By induction, deli Y and linki Y are Laplacian integral, 
and we are done. □ 
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In the case of a shifted simplicial complex A, it is shown in [7] that s ud (A) = <i T (A), where d T 
is the transpose of the vertex-facet degree sequence, and that d T satisfies the simplicial analogue of 
the recursion of Theorem 15.41 

Problem 5.6. Is there an analogous notion of degree sequence for cubical complexes that is related 
to the Laplacian spectrum? 

5.2. Mirroring. Let A be a simplicial complex on vertex set [n]. The mirror of A is the cubical 
complex 

M(A) = {/ e Q n : dir(/) 6 A}. 

Mirroring was first used by Coxeter [6] to study certain generalizations of regular polytopes; see [2j 
section 2.1] for a nice summary of its history and uses. It is not hard to see that the mirroring operator 
takes the class of shifted simplicial complexes to the class of shifted cubical complexes. Mirroring 
also behaves nicely with respect to other basic operations on simplicial and cubical complexes: 

(1) M(Cone(A)) = P(M(A)); 

(2) M(A fe ) = (M(A)) fc ; 

(3) M(dehA) = deliM(A); 

(4) M(linkj A) = link, M (A); 

(5) M(bdA) = bdM(A). 

Here del and link have their usual meanings for simplicial complexes, i.e., 

deli A = {a\i : a G A}, link, A = {a\i : a 6 A, i 6 er}, 

and bdX denotes the union of the boundary faces of every face of X (equivalently, the set of 
non-maximal faces of X). 

Furthermore, mirroring takes near-cones to near-prisms in the following sense. A simplicial near- 
cone A with apex i has the property that bd(deh A) C link, A. In this case, 

bd(del, M (A)) = bd(M (del* A)) = M(bd(deb, A)) C M (linkj A) = link* M(A), 

and so we may construct a near-prism (as in l|26p ) from the cubical complexes deb M(A) and 
linkjM(A). 

Unfortunately, mirroring does not seem to behave nicely with respect to trees or Laplacian eigen- 
values. For instance, mirroring does not preserve the APC property. Thus, the mirror of a simplicial 
complex with spanning trees will not necessarily have a spanning tree. Even the mirror of a pure 
shifted simplicial complex is not necessarily APC. (For example, let A be the graph with vertices 
1,2,3 and edges 12, 13. Then A is contractible, but M(A) is the prism over an empty square - 
a 2-dimensional cell complex that is homotopy equivalent to a circle, hence not APC.) Mirroring 
also does not in general preserve the property of being Laplacian integral. (For example, the pure 
2-dimensional complex A with vertices 1, 2, 3, 4, 5 and facets 124, 125, 134, 135 is a matroid complex, 
hence Laplacian integral by [19] . but M(A) is not Laplacian integral.) 

Although mirroring does not appear useful for enumerating spanning trees, or computing eigen- 
values, it is possible that there are still interesting things to explore, such as the relations between 
eigenvalues of a simplicial complex and its mirror, perhaps especially in special cases such as shifted 
or near-cone/near-prism complexes. 

5.3. Homotopy type. Bjorner and Kalai [3] proved that a simplicial near-cone is always homo- 
topy equivalent to a wedge of spheres, and that the number of spheres in each dimension is easily 
determined from its combinatorial structure. An immediate consequence is that shifted simpli- 
cial complexes are homotopy equivalent to wedges of spheres, with the number of spheres in each 
dimension again easy to describe. 

Given a near- prism X, define 

Bi(X) = deb. X\ link* X. 
By the definition of near-prism, every face in O t Bi(X) or in VBi(X) is a facet of X. 
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Conjecture 5.7. If X is a near-prism in direction i, and deli X is homotopy equivalent to a wedge 
of spheres, 

dekX~\JS r /, 

3 

then X is homotopy equivalent to a wedge of spheres. Specifically, 

.v •- \/ s;! ' . 

3 

where Cij is the number of j- dimensional cells in Bi(X). 

Problem 5.8. Assuming the truth of Conjecture 15.71 give a combinatorial formula for the homotopy 
type of a shifted cubical complex. 

Eran Nevo |21j recently showed us a solution to Problem 15.81 In particular, his combinatorial 
formula confirms Conjecture 15.71 when X is a shifted cubical complex. 

5.4. Extremality. Shifted simplicial complexes derive their name from the existence of shifting 
operators which associate a shifted complex to any simplicial complex. Shifted simplicial complexes 
exhibit certain extremality properties with respect to invariants such as /-vectors, Betti numbers, 
degree sequence and Laplacian eigenvalues: see, e.g., [7J[T7|[I3]. In [3J, Bjorner and Kalai pose the 
problem of developing shifting operators for arbitrary polyhedral complexes. 

Problem 5.9. Is there a natural notion of cubical shifting which associates a shifted cubical complex 
to an arbitrary cubical complex? In what ways are shifted cubical complexes extremal in the class 
of all (proper) cubical complexes? 

6. Duality 

In this section we examine pairs of dual complexes. We now extend the definition of cell complex 
to allow the possibility that the complex contains the empty set as a (— l)-dimensional face. We 
assume the reader is familiar with the basics of matroid theory and refer to [22) for an excellent 
reference. 

Let X and Y be cell complexes. We say that X and Y are dual if there is an inclusion-reversing 
bijection />->■/* from the cells of X to the cells of Y such that, for some d, dim / + dim/* = d for 
every /. (Necessarily, then, dimX = d+ 1 if Y contains the empty set as a face, and dimX = d 
otherwise; similarly, diml" = d+ 1 if X contains the empty set as a face, and diml" = d otherwise.) 
We also require that the boundary maps of X equal the coboundary maps of Y: that is, if we extend 
the duality bijection to all cellular chains by linearity, then 

or, more specifically, 

d Xji : d{X) Ci-i(X) = d*^ d _ l+l : C d ^{Y) -> C d ^ +1 (Y) and 
d* X4 : C t (X) -> CWpO = d Y4 _ i+1 : C d ^{Y) -> C d . i+1 (Y). 

These equalities imply immediately that 

L'x, = L Y , d -i, E'(X; r) = E'^Y; r), (28) 

where • € {ud, du, tot}. 

Proposition 6.1. Suppose that X d and Y d are dual complexes. Let T C Xi and let U = {/* : / G 
Xi\T}. Then the subcomplex Xt — TU^.i is an i-spanning tree of X if and only if the subcomplex 
Yjj — U U Yd—i—i is a (d — i)-spanning tree ofY. 
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Proof. An equivalent statement is the following: the matroicH represented (over any field of charac- 
teristic 0, say Q) by the columns of M — d Xi is dual to the matroid represented by the columns of 
or equivalently by the rows of N — d X i+1 . Since d X i d x i+1 — 0, the column span 
of d Xi (regarded as a vector subspace of QXi) is orthogonal, under the standard inner product, to 
the row span of d Xi+1 . On the other hand, these two subspaces have complementary dimension, 
because 

rankM = |Xi|-dimkerM = |X,|-dimimM = \Xi\-mnkN = |Ai|-rankiV. 
The desired duality now follows from [221 Exercise 2.2.10(h)]. □ 

Duality carries over naturally to the algebraically weighted setting. Specifically, let each cell 
/ € X have an indeterminate weight £/. We define the weighted cellular boundary and coboundary 
maps as in Section 12.61 which give rise to weighted Laplacians and spectrum polynomials as usual. 
Then, we assign weights to the cells of the dual complex Y by the formula 

£/• = I/O- 

It is routine to check that the formulas of (|28p carry over to the weighted setting, that is, 

L x ,i=LY, d -i and E'(X;r)=E' d _ i (Y;r), (29) 

where • <E {ud, du, tot}. 

Furthermore, the matroid represented by dx,i is identical to that represented by d Xi , since we 
have just adjoined the indeterminates {£/ : / G X} to the ground field, and then multiplied the 
rows and columns of the matrix by nonzero scalars (which does not change the matroid structure). 
Therefore, the proof of Proposition l6. II is still valid if we replace the boundary and coboundary maps 
with their algebraically weighted analogues. 

6.1. Spectrum polynomials of matroids. Let M be a matroid on ground set V, and let X be 

the corresponding independent set complex, that is, the simplicial complex on V whose facets are 
the bases of M. Kook, Reiner and Stanton [19j defined the spectrum polynomial of M to be 

Spec M (*;r) = ^«rl< Jl >l, 
i 

where r is the rank function of M; A i— > (A) is its closure operator; and 1\ is a certain subset of / 
defined algorithmically (the details are not needed in the present context). The algorithm depends 
on a choice of a total ordering for V, but Spec M (i;r) does not. By [19l Corollary 18], the spectrum 
polynomial records the Laplacian eigenvalues of X , via the formula 

]T ]T fr x = E tot (X;r,t) = t^r™ Spec M (t; r' 1 ). (30) 
i \es\ ot (x) 

We are going to apply this result to calculate the Laplacian spectra of a complete colorful complex, 
and thereby recover Adin's torsion- weighted count of their simplicial spanning trees [1] . We will see 
that the dual to the n-cube arises as a special case of a complete colorful complex. 

Let oi, . . . , a n be positive integers. The complete colorful complex X = A(ai, . . . , a n ) is defined 
as follows. Let V\ = {i>i,i, . . . , vi, ai }, ...,V n = {v nj i, . . . , v n , an } be n pairwise disjoint vertex sets of 
cardinalities a\, . . . , a n . We will regard each Vi as colored with a different color. Then X is the pure 
(n — l)-dimcnsional simplicial complex on V — V\ U • • • U V n defined by 

X = {fCV: |/nVi|<l Vz}. 

2 When X is a simplicial complex, this matroid is known as a simplicial matroid; see, e.g., \E\. (One could 
analogously define a cellular matroid as a matroid rcpresentable by the columns of some cellular boundary map. This 
is a more general class of matroids; for example, simplicial matroids must be simple, while cellular matroids need not 
be.) 
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Define the color set of an independent set I to be the set C(J) = {i : I n Vi ^ 0}. Also, for a set of 
colors K C [n], abbreviate K = [n]\K and A(K ) := J^i£jf a *- 

Observe that X is a matroid independence complex of a very simple form; the matroid M is just 
the direct sum of n rank-1 matroids whose ground sets are the color classes Vi. In particular, X is 
Cohen-Macaulay [2H pp. 88-89], so H l (X,Z) = for every i < n - 2. 

If n = 2, then X is a complete bipartite graph, while if ai = 2 for every i then X is the boundary 
of an n-dimensional cross-polytope. 

Let W = {171,1, ifyi, ■ ■ • , u n ,i}. For any independent set I C V (i.e., any face of X), the subset Ji 
produced by the algorithm of [TO] is just I\ = I\W; in particular, = A(C{I\)). Therefore, 

Spec M (i;r)=£ £ ftt r l<A>l 

KC[n] IEX: C(I)=K 



= E E 

A'C[n] /£X: C{I)=K 

= E <|K| n( 1 + ( o *- 1 ) r °*)- 



(31) 



A'C[nl fee A 



Applying the Kook-Reiner-Stanton theorem (fSOj) and extracting the coefficient of i s (for — 1 < i < 
n — 1 = dimX), we obtain 



£T(X;r)= E 



KC[n] 
\K\=i+l 



IK' 

\keK 



a k -l) 



e (E n - 1) 1 ^ 

JCA \jGA\J 



AC[n] 
|A-|=i+l 



E (E n c^-D 

jca \jeK\J 



KC[n] 
1X1=1+1 



Therefore, 



n ri A ( j ") 

KC[n] JCA 
|K|=j+l 

Plugging (|33|) into (fTTj) gives 



n 



n-\J\ 
i + \J\ 



JC[n] 
\J\<i+l 



fe-1 



w = n n 



(— l) fc -* +1 (*r — ») 



n-\J\ 
i+l-\J\ 



i=0 JC[n] 
\J\<i+l 



n n b ( j ) 

j=0JC[n] 
\J\=j 



k — 1 

J2 (-i) fe - m (fc-z) 
i=j~i 



n-j 
i + l-j 



(32) 



(33) 



(34) 
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where 



B(J) = A{J) 

The bracketed exponent in (|34l) can be simplified, first by rewriting it in terms of the three quantities 
N = n — j , M = k — j , a = i + 1 — j, and then by some routine calculations (which we omit) using 
Pascal's recurrence: 

fc_1 f n - n \ M /N 

£ (-!)*-<+! ik-i)(l 3 _)= (-1) M (M - a + 1) ' 

i=j-l ^ a=0 

= (-i) M B- i rf' v " lN 

'i\T-2\ _ fn-j-2^ 
M ) \ k-j 

This is precisely the exponent that appears in Adin's formula for Tk(X) [TJ Theorem 1.5]. (Adin's r 
is our ri, and Adin's ci is our j.) In particular, it is zero when j > k. Rewriting (|34[) recovers Adin's 
formula for the tree enumerators of complete colorful complexes: 

Theorem 6.2 (Adin). Let X be the complete colorful complex with color classes of sizes a±, . . . , a n , 
and define B( J) as above. Then 

k 

MX) = n n b{j)^-t). 

j=0 JC[n] 
\J\=J 

In the special case that a, = 2 for all i (so X is the boundary of an n-dimensional cross-polytope), 
the formula (|3ip specializes to 



Spec M (i;r)=^yt fe (l + r 2 ) fc 
and then applying (|30|) . we see that 

E E * <fA = Y,t k ' 1 (f)r 2n - 2k {l + r' 1 ) k = J2 v(. n X 2n - 2 i- 2 (l + r 2 y +1 

i AGsWX„1 k=0 ^ ' j=-l ^ + ' 



i \es\ ot {X n ) k =0 3 

or equivalently 



tot 



The dual complex to an n-dimensional cross-polytope is the n-cube Q n . By formula (|28[) (with 
d = n — 1), we have therefore 

Ef\Q n -r) = i^^r) = Q r 2fc (l + r 2 )- fc , 

giving another proof of (the first formula of) Theorem 13.41 
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